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Summary. — We search for evidence of a light Higgs boson (yl") in tire radiative 
decays of the narrow T(3S') resonance: T(35') — > 7^*^, where A!^ — j- invisible or A^ 

Such an object appears in extensions of the Standard Model, where a light 
CP-odd Higgs boson naturally couples strongly to &-quarks. We find no evidence 
for such processes in a sample of 122 x 10*^ r(35') decays collected by the BABMi 
collaboration at the PEP-II B-factory, and set 90% C.L. upper limits on the product 
of the branching fractions B{T{iS) 7^4°) x invisible) at (0.7-31) x 10"'' 

in the mass range m^o < 7.8 GeV, and on the product B(r(35') ^ 7^*') x ^ 
/i+^i") at (0.25 - 5.2) x 10"^ in the mass range 0.212 < m^o < 9.3 GeV. We also 
set a limit on the dimuon branching fraction of the recently discovered r]i, meson 
B{rib — » jj.'^fj,^) < 0.8% at 90% C.L. The results are preliminary. 

PACS 12.60.Fr - Extensions of electroweak Higgs sector. 
PACS 14.80.Cp - Non-standard-model Higgs bosons. 

PACS 13.66.Fg - Gauge and Higgs boson production in e"e+ interactions. 



INTRODUCTION 



The concept of mass is one of the most intuitive ideas in physics since it is present 
in everyday human experience. Yet the fundamental nature of mass remains one of the 
greatest mysteries in physics. The Higgs mechanism is a theoretically appealing way 
to account for the different masses of elementary particles [1]. The Higgs mechanism 
implies the existence of at least one new particle called the Higgs boson, which is the 
only Standard Model (SM) [2] particle yet to be observed. If it is found, its discovery will 
have a profound effect on our fundamental understanding of matter. A single Standard 
Model Higgs boson is required to be heavy, with the mass constrained by direct searches to 
rriH > 114.4 GeV [3] and rrin 7^ 170 GeV [4], and by precision electroweak measurements 
to ruH = 129t^9 GeV [5]. 

The Standard Model and the simplest electroweak symmetry breaking scenario suffer 
from quadratic divergences in the radiative corrections to the mass parameter of the Higgs 
potential. Several theories beyond the Standard Model that regulate these divergences 
have been proposed. Supersymmetry [6] is one such model; however, in its simplest 
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form (the Minimal Supersymmetric Standard Model, MSSM) questions of parameter 
fine-tuning and "naturalness" of the Higgs mass scale remain. 

Theoretical efforts to solve imattractive features of MSSM often result in models 
that introduce additional Higgs fields, with one of them naturally light. For instance, 
the Next-to-Minimal Supersymmetric Standard Model (NMSSM) [7] introduces a singlet 
Higgs field. A linear combination of this singlet state with a member of the electrowcak 
doublet produces a CP-odd Higgs state A° whose mass is not required to be large. Direct 
searches typically constrain m^o to be below 2mb [8] making it accessible to decays of 
T resonances. An ideal place to search for such CP-odd Higgs would be T ^ 7^", as 
originally proposed by Wilczek [9]. A study of the NMSSM parameter space [10] predicts 
the branching fraction to this final state to be as high as 10~*. 

Other new physics models, motivated by astrophysical observations, predict similar 
light states. One recent example [11] proposes a light axion-like pseudoscalar boson a 
decaying predominantly to leptons and predicts the branching fraction B(T — > 7 a) to be 
between 10^*' 10^"'' [11]. Empirical motivation for a light Higgs search comes from the 
HyperCP experiment [12]. HyperCP observed three anomalous events in the S pfji'^ 
final state, that have been interpreted as a light scalar with mass of 214.3 MeV decaying 
into a pair of muons [13]. The large datasets available at BABMi allow us to place 
stringent constraints on such models. 

If a light scalar ^4*^ exists, the pattern of its decays would depend on its mass. In 
dark matter inspired scenarios, invisible decays could be dominant. For low 

masses m^o < 2mx, relevant for the axion [11] and HyperCP [12] interpretations, the 
dominant decay mode should be ^ — > li^ lJi~ ■ Significantly above the tau-pair threshold, 
^4" T^T^ would dominate, and the hadronic decays may also be significant. 

Preliminary results from search for invisible Higgs decays are described in Ref. [14]. 
This analysis [15] searches for the radiative production of Higgs in Y{3S) decays, which 
subsequently decays into muons: 

T{3S) -/A°; A° ^i+ii- 

The current best limit on the branching fraction B{T 7^*^) with A^ comes 
from a measurement by the CLEO collaboration on T{1S) [16]. The quoted limits on 
B{Y{1S) 7^0) X B{A° H+IJL-) are in the range (1-20) x 10"*^ for m^o < 3.6 GeV. 
There are currently no competitive measurements at the higher-mass T resonances or for 
the values of m^o above the tt threshold. 

In the following, we describe a search for a resonance in the dimuon invariant mass 
distribution for fully reconstructed final state T(3S') ^{ji^ We assume that the 
decay width of the resonance is negligibly small compared to experimental resolution, as 
expected [11, 17] for m^o sufficiently far from the mass of the 7/6 [18]. We also assume 
that the resonance is a scalar (or pseudo-scalar) particle; while significance of any peak 
does not depend on this assumption, the signal efficiency and, therefore, the extracted 
branching fractions are computed for a spin-0 particle. In addition, following the recent 
discovery of the 77^ meson in T(35') decays [18], we look for the leptonic decay of the 
rib through the chain T(3S') ^ 777^, rjb M'^M ■ If the recently discovered state is the 
conventional quark-antiquark rjb meson, its leptonic width is expected to be negligible. 
Thus, setting a limit on the dimuon branching fraction sheds some light on the nature 
of the recently discovered state. We assume T{rib) = 10 MeV, which is expected in most 
theoretical models and is consistent with BABMi results [18]. 
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2. - THE BABAR DETECTOR AND DATASET 

We search for two-body transitions Y{3S) 7^": followed by the decay A° 
invisible [14] or A" /i+^u" [15] in a sample of (121.8 ± 1.2) x 10^ r{3S) decays col- 
lected with the BABAR detector at the PEP-II asymmetric-energy e+e^ collider at the 
Stanford Linear Accelerator Center. The data were collected at the nominal center- 
of-mass (CM) energy Ecm = 10.355 GeV. The CM frame was boosted relative to the 
detector approximately along the detector's magnetic field axis by f3z ~ 0.469. 

We use a sample of 78.5 fb~^ accumulated on T{4S) resonance (T(45') sample) for 
studies of the continuum backgrounds; since T(45) is three orders of magnitude broader 
than T{3S), the branching fraction T(4S') — !■ is expected to be negligible. For 
characterization of the background events and selection optimization we also use a sample 
of 2.4 fb"^ collected 30MeV below the T(3S') resonance. 

The BABAR detector is described in detail elsewhere [19]. We use the GEANT4 [20] 
software to simulate interactions of particles traversing the BABAR detector, taking into 
account the varying detector conditions and beam backgrounds. 

3. - EVENT SELECTION FOR A" ^ ^u+yu" DECAYS 

Wc; select events with exactly two oppositely-charged tracks and a single; c;nergetic 
photon with a CM energy E* > 0.5 GeV. We allow other photons to be present in the 
event as long as their CM energies are below 0.5 GeV. We assign a muon mass hypothesis 
to the two tracks (henceforth referred to as muon candidates) , and require that they form 
a geometric vertex with the Xvtx < 20 (for 1 degree of freedom), displaced transversely 
by at most 2 cm from the nominal location of the e+e" interaction region. We perform 
a kinematic fit to the T{3S) candidate formed from the two muon candidates and the 
energetic photon, constraining the CM energy of the T{3S) candidate, within the beam 
energy spread, to the total beam energy ^/s. We also assume that the T{3S) candidate 
originates from the interaction region. The kinematic fit improves the invariant mass 
resolution of the muon pair. We place a requirement on the kinematic fit Xr(3S) ^ 
(for 6 degrees of freedom), which corresponds to the probability to reject good kinematic 
fits of less than 10^^. The kinematic fit x^j together with a requirement that the total 
mass of the T{3S) candidate is within 2 GeV of ^/s, suppresses background events with 
more than two muons and a photon in the final state, such as cascade decays T(35) — > 
7X6(2P) "f"fT{lS) "f"fii~^fi~ etc. Wc further require that the momentum of the 
dimuon candidate A° and the photon direction are back-to-back in the CM frame to 
within 0.07 radians, and select events in which the cosine of the angle between the muon 
direction and A" direction in the center of mass of is less than 0.88. We reject events 
in which neither muon candidate is positively identified in the muon chamber. 

The kinematic selection described above is highly efficient for signal events. After the 
selection, the backgrounds are dominated by two types of QED processes: "contimium" 
e'^e~ ^jjL^ ii~ events in which a photon is emitted in the initial or final state, and the 
initial-state radiation (ISR) production of the vector mesons J/V', i^i^S), and T{1S), 
which subsequently decay into muon pairs. In order to suppress contributions from 
ISR-produced p° tt+tt" and <j) K~^K~ final states in which a pion or a kaon is 
misidentified as a muon or decays {e.g. through l^-'^'^tJ.), we require that both muons 

are positively identified when we look for A'^ candidates in the range m^o < 1.05 GeV. 
Finally, when selecting candidate events in the rjb mass region m^^ ~ 9.39 GeV, we 
require that no secondary photon above a CM energy of = 0.08 GeV is present in 
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the event; this requirement suppresses decay chains T{3S) 72Xb(2«S') 7172^(16'), in 
which the photon 72 has a typical CM energy of 100 MeV. 

Wc use Monte Carlo samples generated at 20 values of m^o over a broad range 
0.212 < m^o < 9.5 GeV of possible A° masses to measure selection efficiency for the 
signal events. The efficiency varies between 24-44%, depending on the dimuon invariant 
mass. 

4. - EXTRACTION OF SIGNAL YIELDS FOR A" DECAYS 

The invariant mass spectrum for the selected candidates in the T{3S) dataset is shown 
in Fig. 1 (left). We extract the yield of signal events as a function of the assumed mass 
nij^o in the interval 0.212 < m^o < 9.3 GcV by performing a scries of unbinned extended 
maximum likelihood fits to the distribution of tlw^ "reduced mass" 

(1) ruR = ^ml^-Aml . 

The choice of this variable is motivated by the distribution of the continuum background 
from e+e~ 'yn'^iJ,~, which is a smooth function of rriR across the entire range of 
interest, in particular, the region near the kinematic threshold to^^ « 2m^ {ttih w 0). 
Each fit is performed over a small range of mu around the value expected for a particular 
m^o . We use the T(45) sample to determine the probability density functions (PDFs) for 
the continuum background in each fit window, which agree within statistical uncertainties 
with Monte Carlo simulations. We use a threshold (hyperbolic) function to describe the 
background below tur < 0.23 GeV; its parameters are fixed to the values determined 
from the fits to the T{4:S) dataset. Elsewhere the background is well described in each 
limited mu range by a first-order (m^j < 9.3 GeV) or second-order {uir > 9.3 GeV) 
polynomial. 

The signal PDF is described by a sum of two Crystal Ball functions [21] with tail 
parameters on either side of the maximum. The signal PDFs are centered around the 

expected values of rriR = ^Jm^j^a — Arrfi^ and have the typical resolution of 2 — 10 MeV, 
which increases monotonically with m^io. We determine the PDF as a function of m^o 
using a set of high-statistics simulated samples of signal events, and we interpolate PDF 
parameters and signal efficiency values linearly between simulated points. We determine 
the uncertainty in the PDF parameters by comparing the distributions of the simulated 
and reconstructed e+e~ "yiSKJl''l'i J 1^ ^ M"*"/^" events. 

Known resonances, such as J/^/j, i>{^S), and T(IS'), are present in our sample in spe- 
cific intervals of mn, and constitute ■peaking background. We include these contributions 
in the fit where appropriate, and describe the shape of the resonances using the same 
functional form as for the signal, a sum of two Crystal Ball functions, with parameters 
determined from the dedicated MC samples. Wc do not search for A° signal in the 
immediate vicinity of J/ip and tp(2S), ignoring the region of ~ ±40 MeV around J/tp 
(approximately ±5cr) and w ±25 MeV (« ±3cr) around tp{2S). 

For each assunicid value of m^o, we perform a likelihood fit to the mR distribution 
under the following conditions: 

• 0.212 < m,Ao < 0.5 GeV: we use a fixed interval 0.01 < mR < 0.55 GeV. The fits 
are done in 2 MeV steps in m^o. We use a threshold function to describe the 
combinatorial background PDF below rriR < 0.23 GeV, and constrain it to the 
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shape determined from the large T(4S') dataset. For mR > 0.23 GeV, we describe 
the background by a first-order Chebyshev polynomial and float its shape, while 
requiring continuity at mR = 0.23 GeV. Signal and background yields are free 

parameters in the fit. 

• 0.5 < m^o < 1.05 GeV: we use sliding intervals n — 0.2 < niR < ji + 0.1 GeV 
(where ji is the mean of the signal distribution of m^j). We perform fits in 3 MeV 
steps in 7/7 40. First-order polynomial coefficient of the background PDF, signal 
and background yields are free parameters in the fit. 

• 1.05 < m^o < 2.9 GeV: we use sliding intervals — 0.2 < tur < fj, + 0.1 GeV 
and perform fits in 5 MeV steps in m^o. First-order polynomial coefficient of the 
background PDF, signal and background yields are free parameters in the fit. 

• 2.9 < nij^a < 3.055 GeV and 3.135 < 777^0 < 3.395 GeV: we use a fixed interval 
2.7 < iriR < 3.5 GeV; 5 MeV steps in ttj^o. First-order polynomial coefficient of 
the background PDF, signal, J ftp, and background yields are free parameters in 
the fit. 

• 3.4 < m^o < 3.55 GeV: we use sliding intervals /x — 0.2 < tur < jj, + 0.1 GeV 

and perform fits in 5 MeV steps in 771^0 . First-order polynomial coefficient of the 
background PDF, signal and background yields are free parameters in the fit. 

• 3.55 < 777^0 < 3.66 GeV and 3.71 < 777^0 < 4.0 GeV: we use fixed interval 3.35 < 
niR < 4.1 GeV; 5 MeV steps in m^o. First-order polynomial coefficient of the 
background PDF, signal, il){2S), and background yields are free parameters in the 
fit. 

• 4.0 < 777 AO < 9.3 GeV: we use sliding intervals /u - 0.2 < niR < -t- 0. 1 GeV; 5 MeV 
steps in niAo- First-order polynomial coefficient of the background PDF, signal 
and background yields are free parameters in the fit. 

• 7/5 region (ttj^^ = 9.390 GcV): wc use a fixed interval 9.2 < rriR < 9.6 GeV. We 
constrain the contribution from e+e~ —>■ 7isr^(1<S') to the expectation from the 
T{AS) dataset (436 ±50 events). Background PDF shape (second-order Chebyshev 
polynomial), yields of T{3S) jXb(2P) j-fT{lS), signal Y{3S) 7% events, 
and background yields are free parameters in the fit. 

The step sizes in each interval correspond approximately to the resolution in tuao ■ 

5. - SYSTEMATIC UNCERTAINTIES for n+ij,- 

The largest systematic uncertainty in B{T{3S) — > jA^) comes from the measurement 

of the selection efficiency. Wc compare the overall selection efficiency between the data 
and the Monte Carlo simulation by measuring the absolute cross section da/druR for 
the radiative QED process e+e~ — > 7/i+/x~ over the broad kinematic range < tur < 
9.6 GeV, using a sample of 2.4 fb"^ collected 30 MeV below the T{3S). We use the 
ratio of measured to expected cross sections to correct the signal selection efficiency as a 
function of 777^0 . This correction reaches up to 20% at low values of rn^o . We use half 
of the applied correction, or its statistical uncertainty of 2%, whichever is larger, as the 
systematic uncertainty on the signal efficiency. This uncertainty accounts for effects of 
selection efficiency, reconstruction efficiency (for both charged tracks and the photon), 
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trigger efficiency, and the uncertainty in estimating the integrated luminosity. We find 
the largest difference between the data and Monte Carlo simulation in modeling of muon 
identification cflicicncy. 

We determine the uncertainty in the signal and peaking background PDFs by com- 
paring the data and simulated distributions of e+e~ 7isR"^/V' events. We correct for 
the observed 24% difference (5.3 McV in the simulations versus 6.6 MeV in the data) 
in the width of the mR distribution for these events, and use half of the correction to 
estimate the systematic uncertainty on the signal yield. This is the dominant systematic 
uncertainty on the signal yield for m^o > 0.4 GeV. Likewise, we find that changes in 
the tail parameters of the Crystal Ball PDF describing the J/tp peak lead to variations 
in event yield of less than 1%. We use this estimate as a systematic error in the signal 
yield due to uncertainty in tail parameters. 

We find excellent agreement in the shape of the continuum background distributions 
for ruR < 0.23 GeV between r{3S) and r(45) data. We determine the PDF in the 
fits to T{4S) data, and propagate their uncertainties to the T{3S) data, where these 
contributions do not exceed (t{B) = 0.3 x 10^^. For the higher masses mn > 0.23 GeV, 
the background PDF parameters are floated in the likelihood flt. 

We test for possible bias in the fitted value of the signal yield with a large ensemble 
of pseudo-experiments. For each experiment, we generate a sample of background events 
according to the number and the PDF observed in the data, and add a pre-determined 
number of signal events from fully-reconstructed signal Monte Carlo samples. The bias is 
consistent with zero for all values of m^o , and we assign a branching fraction uncertainty 
of a{B) = 0.02 X 10~^ at all values of m^o to cover the statistical variations in the results 
of the test. 

The uncertainties in PDF parameters of both signal and background and the bias 
uncertainty affect the signal yield (and therefore significance of any peak); signal efficiency 
uncertainty does not. The effect of the systematic uncertainties on the signal yield is 
generally small. The statistical and systematic uncertainties on the branching fraction 
B{T{3S) 7^*^) as a function of m^o are shown in Fig. 1 (right). 

6. - RESULTS for ^l+ ^l- 

For a small number of fits in the scan over the T{2>S) dataset, we observe local 
likelihood ratio values S of about 3fT. The most significant peak is at rn^o = 4.940 ± 
0.003 GeV (likelihood ratio value S = 3.0, including systcmatics; B = (1.9 ± 0.7± 0.1) x 
10~^). The second most-significant peak is at m^o = 0.426 ± 0.001 GeV (likelihood 
ratio value S = 2.9, including systematics; B = (3.1 ± 1.1 ± 0.3) x 10~^). The peak 
at 171^0 = 4.940 GeV is theoretically disfavored (since it is significantly above the r 
threshold), while the peak at m^o = 0.426 GeV is in the range predicted by the axion 
model [11]. However, since our scans have 0(2000) rn^o points, we should expect several 
statistical fluctuations at the level of 5 3, even for a null signal hypothesis. At least 
80% of our pseudo-experiments contain a fluctuation with iS = 3(7 or more. Taking this 
into account, we conclude that neither of the above-mentioned peaks are significant. 

Since we do not observe a significant excess of events above the background in the 
range 0.212 < m^o < 9.3 GeV, we set upper fimits on the branching fraction B{T{3S) 
jA°) X B{A^ /i+/U~). We add statistical and systematic uncertainties (which include 
the additive errors on the signal yield and multiplicative uncertainties on the signal 
efficiency and the number of recorded T{3S) decays) in quadrature. The 90% C.L. 
Bayesian upper limits, computed with a uniform prior and assuming a Gaussian likelihood 
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Fig. 1. - Distribution of the dimuon invariant mass m^+^- in the T(3S') data is shown on the 
left. Statistical and systematic uncertainty on the product of branching fractions B{T{?>S) —> 
yA°) X B{A'^ fj,'^fj,~) are shown on the right as a function of m^o, extracted from the fits to 
the T{3S) data. Statistical errors are shown as red dot-dashed line, systematic uncertainties are 
shown as blue dotted line, and the total uncertainty, computed as a quadrature sum of statistical 
and systematic errors, is the solid black line. The shaded areas show the regions around the 
J/ip and ip(2S) resonances excluded from the search. 



function, are shown in Fig. 2 (right), as a function of mass m^o. The limits fluctuate 
depending on the central value of the signal yield returned by a particular fit, and range 
from 0.25 x 10^*^ to 5.2 x 10"^ 

We do not observe any significant signal at the HyperCP mass, m^o = 0.214 GeV. 
We find S(T(3S') 7^10(214)) = [O.Utafi ± 0.17) x 10-^ and set an upper limit 
S(T(3S') 7^0(214)) < 0.8 x 10-^ at 90% C.L. 

From a fit to the ryf, region, we find B[T[iS) — > 777^) x B{rj}j — > M^/^^) = (0.2 ± 
3.0 ± 0.9) X 10~^, consistent with zero. Taking into account the BABAR measurement of 
K(T(3S') ^ 7r?h) = (4.8 ± 0.5 ± 1.2) x 10^^ we can derive B{r]i, -> ) = (0.0 ± 0.6 ± 

0.2)%, or an upper limit /B(?/h ii^ji^) < 0.8% at 90% C.L. This is consistent with 
expectations from the quark model. All results above are preliminary. 

The limits we set [15] are more stringent than those recently reported by the CLEO 
collaboration [16]. Our limits rule out much of the parameter space allowed by the light 
Higgs [10] and axion [11] models. 

7. CONCLUSIONS 

We find no evidence for light Higgs boson in a sample of 122 x 10^ T{3S) decays 
collected by the BABAR collaboration at the PEP-II B-factory, and set 90% C.L. upper 
limits on the product of the branching fractions B{T{3S) 7^*^) x B{A'^ invisible) 
at (0.7 — 31) X 10^^ in the mass range to^o < 7.8 GeV [14] and on the product 
B(T{SS) 7AO) X B{A° fJ.+^J-~) at (0.25 - 5.2) x 10"^ in the mass range 0.212 < 
m^o < 9.3 GeV [15]. We also set a limit on the dimuon branching fraction of the re- 
cently discovered rjh meson B{?]b M^M ) < 0.8% at 90% C.L. [15]. The results are 
preliminary. 
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Fig. 2. - Distribution of the likelihood ratio variable 5 with additive systematic uncertainties 
included for the fits to the T(35) dataset, overlayed with a blue curve showing the Gaussian fit 
with fixed fi = and a = 1, is shown on the left. Upper limits on the product of branching 
fractions B{T{3S) 7^") x B{A° — > fj,'^ fj, ) as a function of m^o from the fits to T{SS) data are 
shown on the right. The shaded areas show the regions around the J/tp and '4>{2S) resonances 
excluded from the search. 
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